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Abstract 

Future generation of interferometric gravitational wave detectors is hoped to provide 
accurate measurements of the final stages of binary inspirals. The sources probed 
by such experiments are of extragalactic origin and the observed chirp mass is the 
intrinsic chirp mass multiplied by (1 + z) where z is the redshift of the source. 
Moreover the luminosity distance is a direct observable is such experiments. This 
creates the possibility to establish a new kind of cosmological tests, supplementary 
to more standard ones. 

Recent observations of distant type la supernovae light-curves suggest that the 
expansion of the universe has recently begun to accelerate. A popular explanation of 
present accelerating expansion of the universe is to assume that some part of the 
matter-energy density is in the form of dark component called "the quintessence" 
with the equation of state pQ = wpq with w > —1. In this paper we consider 
the predictions concerning observations of binary inspirals in future LIGO type 
interferometric experiments assuming a "quintessence cosmology". In particular 
we compute the expected redshift distributions of observed events in the a priori 
admissible range of parameters describing the equation of state for the quintessence. 
We find that this distribution has a robust dependence on the cosmic equation of 
state. 
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1 Introduction 



Recent distance measurements from high-redshift type la supernovae ^ suggest that 
the universe is presently accelerating its expansion. A popular explanation of this phe- 
nomenon is to assume that considerable amount Qq ~ 70% of the matter-energy density 
is in the form of dark component called "the quintessence" characterised by the equation 
of state pq = wpq with w > — 1 ^, The evidence for spatially flat universe, rein- 
forced by recent cosmic microwave background experiments BOOMERANG and MAX- 
IMA 0] calls for an extra unclustered dark component. Within the standard cold dark 
matter (CDM) scenario only about 0.2 < VLcdm < 0.4 can be clustered in order to be in 
agreement with galactic rotation curves, abundance of galaxy clusters, gravitational lens- 
ing or large scale velocity fields. Moreover the accelerated expansion of the universe can 
be achieved with extreme forms of matter. Hence this extra component (quintessence) 
should be similar to the cosmological constant but is allowed to have its own temporal 
dynamics. Many current models of dark matter in general and of quintessence in partic- 
ular [§, invoke the concepts from particle physics. Particle physics, however, gives little 
guidance as to concrete models of quintessence. Therefore it has been proposed in p that 
future supernova surveys may allow reconstructing the quintessential equation of state. 
In this paper we shall contemplate the feasibility of constraining the cosmic equation of 
state from the gravitational wave experiments in a similar vein as proposed in ||^. 

Laser interferometric gravitational wave detectors developed under the projects LIGO, 
VIRGO and GEO600 are expected to perform a successful direct detection of the gravi- 
tational waves. Inspiralling neutron star (NS-NS) binaries are among the most promising 
astrophysical sources for this class of experiments |T^. Besides quite obvious benefits from 
seeing gravitational waves "in flesh" and providing valuable information about dynamical 
processes leading to their generation inspiralling binaries have one remarkable feature. 
Namely, the luminosity distance to a merging binary is a direct observable quantity easy 
to obtain from the waveforms. This circumstance made it possible to contemplate a possi- 
bility of accurate measurements of cosmological parameters such like the Hubble constant. 



or deceleration parameter |]TT], |12| |13[. In particular it was pointed out by Chernoff and 
Finn |]TT| how the catalogues of inspiral events can be utilised to make statistical infer- 
ences about the Universe. In the similar spirit we discuss in this paper the possibility to 
constrain the quintessence equation of state from the statistics of inspiral gravitational 
wave events. 



2 Cosmological model 

We shall consider a class of fiat quintessential cosmological models. The spatially fiat 
Universe has recently received a considerable observational support from the measure- 
ments of the position of first acoustic peak at / ~ 200 in baloon experiments BOOMERANG 
and MAXIMA This class is parametrized by two quantities: Qq and Qq, where 

Qo = p/ Per = 2 ° denotes the current matter density as a fraction of critical den- 

3-Ho 

sity for closing the Universe, VLq is analogous fraction of critical density contained in 
the quintessence and these two sum up to the value one. The equation of state for the 
quintessence is assumed in a standard form: p = wp where w > —1. This form of the 
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equation of state is very general in the sense that it contains the well known constituents 
of the universe as special subclasses. For example w = —1 corresponds to the cosmologi- 
cal constant A, w = - the dust matter, w = —1/3 - cosmic strings and w = —2/3 the 
domain walls. 

Non Euclidean character of the space-time is reflected in distance measures. For the 
introduction to observational cosmology and the problems of distances in non-euclidean 
spaces (see e.g. [0). In order to fix the notation for further use, let us introduce an 
auxihary quantity T>{z): 



Viz) = ^Qoil + zy + Qq{1 + (1) 

As it is well known, one can distinguish three types of distances: 
(i) proper distance: 

Hq Jo V[w) h Jo V[w) h 



(ii) angular distance: 



dA{z) = —^dMiz) = -^^dMiz) (3) 
1 + z 1 + z h 

(iii) luminosity distance: 

dLiz) = (1 + z)dMiz) = (1 + z)^dM{z) (4) 

h 

As usually z denotes the redshift, h denotes the dimensionless Hubble constant i.e. 
Hq = h X 100 km/sMpc and d^ = S. x Gpc is the Hubble distance (radius of the 
Hubble horizon). The quantities with an overbar have been defined by factoring out the 
dependence on the Hubble constant from respective quantities. In the further discussion 
we will explore the following models: 

(fio, 1)q) = {(0.2,0.8); (0.3,0.7); (0.4,0.6)} 

with the w coefficient equal to w = {0, —0.2, —0.4, —0.6, —0.8, —1.} 

/^From the observational point of view in the light of constraints from large scale 

structure and cosmic microwave background anisotropics, the 95% confidence interval 

estimates give 0.6 < < 0.7 and -1. < < -0.6 ^ [T^. 

However we retain the full spectrum of a priori possible quintessential equations of state 

in order to illustrate the discriminating power of the gravitational wave data discussed in 

this paper. 



3 Redshift distribution of observed events 

The gravity wave detector would register only those inspiral events for which the signal- 
to-noise ratio exceeded certain threshold value po Ijn], which is estimated as po = 8- 
for LIGO-type detectors. An intrinsic chirp mass Aio = p^^^M^/^, where p and M denote 
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the reduced and total mass, is the crucial observable quantity describing the inspiralling 
binary system. The observed chirp mass Ai{z) = (1 + z)Aio scales with the redshift and 
therefore can be used to determine the redshift to the source (there is strong evidence 
that the mass distribution of neutron stars in binary systems is sharply peaked around the 
value 1.4 Mq). Because the luminosity distance of a merging binary is a direct observable 
easily read off from the waveforms one has a possibility to determine the precise distance - 
redshift relation and hence to estimate the Hubble constant pDI, ITBl . For a given detector 



and a source the signal-to-noise ratio reads 




P{Z) = 86-^ — ^ Cifmax), (5) 



where tq is a charateristic distance scale, depending on detector's sensitivity, tq ~ 355 Mpc 
for advanced LIGO detectors, is the luminosity distance to the source, Cifmax) is a 
dimensionless function describing the overlap of the signal with detector's bandwidth. 

The adiabatic inspiral signal terminates when the binary system reaches the innermost 
circular orbit (ICO). The corresponding orbital frequency is fico and fmax corresponds 
to observed (i.e. redshifted) fico 

fico no Hz f 2.8 Mq 



f"^-- l + z l + z [ M 
so fmax ~ 710 Hz foT ucutrou star binaries. It is argued that Cifmax) ~ 1 for LIGO/VIRGO 



interferometers 11, 29 . 



Let us denote by no the local binary coalescing rate per unit comoving volume. One 
can use "the best guess" for local rate density tiq ^ 9.9 h 10~^ Mpc~^yr~^ as inferred 
from the three observed binary pulsar systems that will coalesce in less than a Hubble 
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Source evolution over sample is usually parametrized by multiplying the coalescence rate 
by a factor ri{z) = (1 -|- z)^, i.e. h = uq [1 + z)"^ ri{z) where the (1 -|- z)^ factor accounts 
for the shrinking of volume with z and the time dilation of burst rate per unit time. The 
cosmological origin of gamma-ray bursts (GRBs) has been confirmed since discoveries 
of optical counterpart of GRB 970228 |^ and the measured emission-line redshift of 
z = 0.853 in GRB 970508 ^T|. It has also been known for quite a long time that 
cosmological time dilation effects in BATSE catalogue suggest that the dimmest sources 
should be located at 2; ~ 2 [E^. Consequently several authors tackled the question of 



source evolution in the context of gamma-ray bursts. Early estimates of p3[ and Piran 
indicated that BATSE data could accommodate quite a large range of source density 
evoultion (from moderate negative to positive one). Later considerations by Horack et 
al. indicated that if 2; = 2 is indeed the limiting redshift then a source population 
with a comoving rate density n{z) ~ (1 -|- z)^ with 1.5 < (3 < 2 is compatible with 
BATSE data. Later on Totani ||26| considered the source evolution effects and based his 
calculations on the realistic models of the cosmic star formation history in the context of 
NS-NS binary mergers. Comparison of the results with BATSE brightness distribution 
revealed that the NS-NS merger scenario of GRBs naturally leads to the rate evolution 
with 2 < j3 < 2.5. We shall therefore take the source evolution effetcs into account in 
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our further considerations. One should stress, however that NS-NS merger scenario is by 
no means the unique explanation of gamma-ray bursts. Recently the so called collapsar 
model became popular |^. The idea that at least some of gamma-ray bursts are related 
to the deaths of massive stars is supported by the observations of afterglows in GRB 
970228 and GRB 980326 [?]. Therefore we will not prefer any specific value of evolution 
exponent D but instead we will try to illustrate how strongly and in which direction does 
the source evolution affect our ability to discriminate between different quintessential 
equations of state. 

The relative orientation of the binary with respect to the detector is described by the 
factor B. This complex quantity cannot be measured nor assumed a priori. However, its 



probability density averaged over binaries and orientations has been calculated [^] and 
is given by a simple formula: 



Pern 



56(4 
0, 



- 9)7256, 
otherwise 



if < e < 4 



(7) 



The rate - — ^'^^ at which we observe the inspiral events that originate in the 



dz 

redshift interval [z, z + dz] is given by 
dN{> po) no 



d 



dz 



r]{z) Airdlj^dMiz) Ce{x) 
1 + z dz 



An 



'd 



H 



ho dl[{z) 



h 



1 + z Viz) 



CM 



(8) 



where Cq{x) = Pe{Q)dQ denotes the probability that given detector registers inspiral 
event at redshift Zs with p > po- The quantity Ce{x) can be calculated as 



Ce(x) = (l + x)(4-x)V256 
for X > 4 



for < X < 4 



(9) 



where 31 : 



X 



and 



A := 0.4733 



hA 



'l + z 



,7/6 



dA{.z) 
du/h 



— {l + zf"dM{z) 




Mo ^ 
1.2 Mr,, 



5/6 



(10) 



(11) 



^po/ V355 Mpc^ 

Figure 1 shows the expected detection rate of inspiralling events for the cosmological 
model with (f2o = 0.3 , VLq = 0.7) assuming no source evolution and covering the full range 
of a priori possible quintessential equations of states. It has been obtained by numerical 
integration of the formula (^. The predictions for other realistic proportions of flo and Qq 
are almost indistinguishable at the level of detection rates, so the Fig.l is representative for 
the whole class of models considered. The effect of source evolution on the detection rate 
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is summarised in Fig. 2. For transparency only one member (corresponding to Wq = —0.8) 
of each family of curves (as in Fig.l) is shown for different values of the evolution exponent 
D. 

The method of extracting the cosmological parameters advocated by Finn and Cher- 
noff 1^ makes use of the redshift distribution of observed events in a catalogue composed 



of observations with the signal-to- noise ratio greater than the threshold value po- There- 
fore it is important to find this distribution function for different quintessence models. 
The formula for the expected distribution of observed events in the source redshift can 
be easily obtained from the equation (||): 




Vi^) Ce{x) (12) 



^(> Po) V ^ / l + z ' V{z 



The summary of numerical computations for the cosmological quintessence models 
considered based on the formulae (|12|) and (|^) are given in figures Fig. 3 and Fig.4. Fig.3 
illustrates the P(z, > po) distribution function for the (fio = 0.3 , VLq = 0.7) cosmological 
model with different quintessential equations of state. For the purpose of obtaining the 
Figures 3 and 4 we have assumed the dimensionless Hubble constant equal to h = 0.65 



as suggested by independent observational evidence (e.g. SNe la in HST project |^ 
or multiple image quasar systems [^). On Fig.4 the distribution functions for different 
cosmological models with the quintessence field with w = —0.8 have been plotted together. 
Fig. 5 shows the distribution functions for different evolutionary exponent in the {Qq = 
0.3 , fig = 0.7) model with w = -0.8. 



4 Results and discussion 

It is clear from Figure 1 that different quintessential cosmologies (singled out by w pa- 
rameter in the equation of state) give different predictions for annual inspiral event rate 
to be observed by future interferometric experiments. Unfortunately , this difference is 
too small to be of observational importance. Moreover, as already pointed out, there 
exists a degeneracy in terms of cosmological models (labelled by the value of Qq and Qq). 
There is however a difference between detection rates corresponding to different values of 
evolutionary exponents as displayed in Figure 2. 

Figure 3 shows that there is a noticeable difference in predicted event redshift dis- 
tribution functions P{z, > po) for different values of the cosmic equation of state within 
given cosmological model (labelled by the values of Qq and fig). The spread between 
different cosmological models for a given quintessence equation of state is much smaller 
as seen from the Fig.4. This is a reflection of above mentioned effective degeneracy with 
respect to values of Q parameters. Hopefully this degeneracy can be broken by inde- 
pendent estimates of and Qq parameters in other studies (cluster baryons estimates, 
Lya forest surveys, large scale structure or CMBR). The spread of redshift distribution 
functions attributed to evolutionary effects is smaller as shown in Fig. 5 and has a slightly 
different character - the distribution function is shifted toward increasing redshifts when 
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the evolutionary exponent changes from positive to negative value. This may to some 
extent mimic the effect of cosmic equation of state, but it should in principle be possible 
to disentangle - at least to a certain degree from the complementary information about 
the detection rates. As can be seen from Fig.2 the magnitudes of observed event rates 
for different evolutionary exponents D are clearly distinct, at least for the range of the 
Hubble constant suggested by independent cosmological evidence [|^. 

The redshift distribution P(z, > po) is in fact inferred from observed chirp mass distri- 
bution. Therefore it can in principle be distorted by the intrinsic chirp mass distribution. 
Theoretical studies of the neutron star formation suggest that masses of nascent neutron 



stars do not vary much with either mass or composition of the progenitor Also the 



mass estimates of observed binary pulsars suggest that there are good reasons to assume 
a negligible spread of intrinsic chirp mass (as it was done in the present paper). Moreover 
any intrinsic distribution of mass would be expected as symmetric, whereas the redshift 
distribution (of cosmological origin) has certain amount of asymmetry. 

In conclusion one can hope that the catalogues of inspiral events gathered in future 
gravitational waves experiments can provide helpful information about the quintessence 
equation of state complementary to that obtained by other techniques. Even though the 
most straightforward way of making inference about cosmic equation of state would come 
from future supernovae surveys it would be good to have in mind alternative ways of 
reaching the same goal such as the one proposed in the present paper. 
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Figure Captions 



Figure 1 

The detection rate prediction for the advanced gravity wave detectors i.e. with 

signal-to-noise threshold po — 8- ^-nd probing distance tq = 355 Mpc. corresponding 
to quintessence cosmology with different equations of state. 



Figure 2 

The detection rate prediction for the advanced gravity wave detectors i.e. with signal- 
to-noise threshold po = 8. and probing distance ro = 355 Mpc. corresponding to Qq — 
0.3, Qq = 0.7 quintessence cosmology with Wg — —0.8 for different values of evolutionary 
exponent D. 



Figure 3 

Redshift distribution of observed events in the cosmological model with flo = 0.3, flq ~ 
0.7 for different quintessential equations of state. 



Figure 4 

Redshift distribution of observed events in the cosmological quintessence model with 
w — —0.8 Different cosmological models have been plotted collectively. 



Figure 5 

Redshift distribution of observed events in the cosmological model with Qq = 0.3, Q,q — 
0.7 with Wq — —0.8 quintessence for different values of evolutionary exponents D. 
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